Quantum transport with coupled cavities on the Apollonian network 
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We study the dynamics of atomic and photonic excitations in the Jaynes-Cummings-Hubbard 
(JCH) model when each cavity occupies the nodes of an Apollonian network (AN). Such structure 
shows small-world and scale-free properties, and has been investigated within many physical models. 
By numerically diagonalizing the system Hamiltonian in the single-excitation subspace, we evaluate 
the time evolution of an initial state prepared as an even superposition of both atomic and pho- 
tonic modes fully localized at a specific node. For the large hopping regime we provide a detailed 
description of the quantum transport properties of the AN and show that the complex interplay 
between the many cavity-cavity degrees of freedom induces both extended and localized states vary- 
ing periodically. The excitation is most likely to be found in the initial node and its propagation 
strongly depends on the initial conditions. We also discuss the strong atom-cavity coupling regime, 
where atomic and photonic modes propagates identically, and the JCH regime, where the system is 
allowed to roam between atomic and photonic degrees of freedom, since the hopping parameter and 
the atom-cavity coupling strength are of the same order. In this regime, different cavities contribute 
either to the atomic or photonic component mostly, depending on the initial conditions. 

PACS numbers: 42.50.Pq, 42.50.Ex, 03.67.-a, 89.75.Da 



I. INTRODUCTION 

Cavity quantum electrodynamics systems have been 
widely considered as a suitable choice for implementa- 
tion of quantum information processing schemes [H-dj]. 
A very simple cavity system consists of a two-state atom 
coupled to a single mode of an electromagnetic field in- 
side an ideal cavity. Such system is described by the well 
known Jaynes-Cummings (JC) model which is ex- 
actly solvable and has been useful for investigating fun- 
damental aspects of the interaction between light and 
matter. Further generalizations of the JC model include 
multiple atoms, AT- level atom, dissipation, and so forth 

Experimental advances in optical microcavities |6|, 
photonic crystals Q , and superconducting devices Q has 
brought interest to Hubbard-like models where photons 
can hop through an array of coupled-cavity systems. In 
articular, the Jaynes-Cummings-Hubbard (JCH) model 
[Tol | provides a new framework for studying many-body 
phenomena such as phase transitions where strongly cor- 
related photons play the role. Most of recent work has 
focused on the Mott insulator-superfluid quantum phase 
transition of polaritons (combined atomic and photonic 
excitations) |ll| - [l3j . 

Coupled-cavity systems allow the control and measure- 
ment of individual nodes. Therefore, it also provides a 
reliable setup for distributed quantum information pro- 
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cessing [l4| , which requires entanglement generation 
and quantum state transfer between distant nodes in a 
network [l6j]. Recently, Ogden et al. [13] studied quan- 
tum state transfer between two coupled cavities, each 
one containing a single two-level atom. In terms of dclo- 
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calized field modes, they have derived the dynamics on 
the single-excitation subspace for a variety of interaction 
regimes focusing on the atomic state transfer and found 
that, by setting the right parameters, one can achieve a 
high-fidelity state transfer without populating the field 
mode, thus avoiding photon loss. Makin et al. [l8| con- 
sidered a one-dimensional array of cavities and dealing 
with both atomic and photonic excitation dynamics they 
discussed the localized and delocalized behaviour of the 
system. They have also shown that by parabolically cou- 
pling the cavities, or encoding the initial state as a gaus- 
sian wave-packet, one can avoid pulse dispersion. In [l9| . 
Dong et al. discussed the binary transmission dynamics 
in a chain of coupled cavities, each one containing a tree- 
level atom and showed that a certain class of coding and 
decoding protocol can improve the transmission fidelity. 

These results indicate that even in the single-excitation 
subspace the JCH model presents a rich dynamics due 
to the interplay between atomic and photonic degrees 
of freedom. Thus it would be interesting to apply this 
model to more complex structures. Lately, a class of 
networks which are neither completely regular nor com- 
pletely random, named small-world networks [20| , where 
the average length of the shortest path I between two 
nodes increases logarithmically with the network size N 
(I cc In AT), has been widely studied within many fields, 
from social and biological systems [20|. [2ll to classical 
and quantum transport dynamics In particular, 

a class of complex networks called Apollonian networks 
(AN) (25|, which are simultaneously small- world, scale- 
free (that displays a power-law degree distribution), and 
can be embedded in Euclidean space, has brought some 
attention. The free-electron gas (26[, magnetic models 
27j , tight-biding models [28] , correlated electron systems 
291 ] , and quantum walks [30] in such networks have been 
investigated. AN's arise from the problem of space-filling 
packing of spheres (3lj . It can be deterministic or ran- 
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FIG. 1. The second generation (n = 2) Apollonian packing 
(solid lines) and the corresponding AN (dotted lines). 



dom, and depends on the initial configuration. The most 
common packing starts with three mutually touching cir- 
cles. The first generation, n — 1, is achieved by insert- 
ing a maximal circle in the interstice bounded by the 
three initial circles. Further generations are obtained 
by repeating this process for every empty space. This 
packing can be mapped to the AN as shown in Fig. [T] 
At a given generation, the number of nodes is given by 
N(n) = (3 n + 5)/2, and the number of connections by 
C{n) = (3" +1 +3)/2 (n = 0,1,2,...). In each generation a 
new group of nodes sharing the same degree 7 is created. 
More precisely, there are 3"" 1 , 3™~ 2 , 3™~ 3 , 3 2 , 3, 1, and 
3 nodes with degree 7 = 3, 3 • 2, 3 • 2 2 , 3 • 2™" 2 , 3 • 2 n ~\ 
and 2" + l, respectively, where the last number represents 
the three corner nodes. The central node, often referred 
as the hub node, has the largest degree. 

In this paper, we study the dynamics of the JCH 
model defined on the AN. We numerically diagonalizc 
the Hamiltonian in the single-excitation subspace and 
and solve Schrodinger equation thus obtaining the time 
evolution of atomic and photonic occupation probabil- 
ity distribution for fully localized initial states prepared 
in an even superposition of atomic and photonic modes. 
We consider three different energy regimes. For large 
inter-cavity hopping, the photonic component propagates 
practically free, without significantly exchange energy 
with the atomic modes in short time scales. In this 
regime we analyse the photonic component transport 
through the AN when the initial state is prepared in dif- 
ferent nodes and show how it changes by varying the net- 
work size. When the atom-cavity coupling strength dom- 
inates, we have a similar propagation dynamics but with 
both components travelling simultaneously. For compa- 
rable hopping and atom-cavity coupling, the dynamics 
turns out to be more complex since the system can roam 
between atomic and photonic degrees of freedom, thus 
leading to different dynamics. 



In what follows, Sec. HH we introduce the JCH model. 
In Sec. IIIII we present our numerical results for the time 
evolution at different energy regimes. Finally, in Sec. IIV| 
we draw our conclusions. 



II. THE JAYNES-CUMMINGS-HUBBARD 
MODEL 

Let us consider a system of coupled cavities where each 
one occupies the nodes of an Apollonian network. A two- 
level system {\g) , |e)} coupled to a single field mode can 
described by the JC Hamiltonian [3] (in the rotating wave 
approximation) , 

h JC = W-^ + W/a t^ + p( a + a . + fr - a t) ! (1 ) 

where a z \g) = - \g) and u z \e) = |e), cr+ {aj) and a] (a ) 
are, respectively, the atomic and photonic raising (lower- 
ing) operators for cavity j, /3 is the atom-cavity coupling 
strength, to a the atomic transition frequency, ujf the field 
mode frequency, and we set h = 1 for convenience. If 
we include photon evanescent hopping between nearest- 
neighbour cavities (tight-binding approximation) , we can 
represent the whole system by the JCH Hamiltonian 

N(n) N(n) 

# = £ />f - £ ^^i, (2) 
3=1 »>j=i 

where k is the cavity-cavity coupling strength, N(n) 
number of nodes for a fiven AN generation n, and Ay 
the adjacency matrix elements. We now restrict the sys- 
tem basis to the single-excitation subspace. These exci- 
tations can be photonic or atomic, that is \Q) \g, 1) or 
|Q)|e,0), respectively, where Q £ {1, N(n)} speci- 
fies the cavity. Then the whole system state is written 
as the tensor product of every single-cavity state, where 
one excitation must be conserved (the other non-excited 
cavities are in the |g,0) state), so the Hilbcrt space has 
D = 2N(n) = 3" + 5 dimensions. In such basis, the ma- 
trix form of the Hamiltonian ([2]) can be simply expressed 
by 

H lcxc = y Jjv(n) ® Z + pI N (n) ®X- kA„ ® — ^ - , (3) 

where A = ujf — uj a is the detuning between photonic and 
atomic modes, I m is the m x m identity matrix, X and 
Z are the usual Pauli matrices, and A„ the adjacency 
matrix. The ground state ®q£" l ' ) \Q) \g, 0) is not included 
in this subspace, that is, the system is already set with 
the single-excitation. 

The AN does not take geometrical deformations into 
account so we are only interested in the connections be- 
tween the nodes. By considering uniform cavity-cavity 
coupling, the adjacency matrix A„ for n = 2, for exam- 



3 




modes (i.e., a single-cavity eigenstate for A = yj), 
|^fc(0)) = |fc) ® (\g, 1) + |e, 0))/y/2 at node (cavity) k, we 
write 



(5) 



and 



(6) 



as the probability of finding both photonic and atomic 
excitations, respectively, at node I in time t, where we 
have denoted \£,g) = \l) \g, 1) and \t,e) = \f) |e, 0). 



FIG. 2. (Color online) The Apollonian network for n = 4 
(43 nodes). The first generation is made up by nodes 1 to 
4, the second one includes 5 to 7, the third, 8 to 16, and the 
fourth, 17 to 43. The corner nodes were drawn bigger for a 
clearer visualization. This is the node labelling we will adopt 
through the paper. 
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(see Fig. [T]). In order to discuss our results in next sec- 
tion, we need the appropriate identification of each node 
(cavity), then Fig. [5] provides the node numbering which 
we are going to consider herein, up to the 4th AN gener- 
ation. 



III. NUMERICAL RESULTS 

In this section we display the results for time-evolution 
of the JHC model in the AN obtained by numerical diag- 
onalization of the Hamiltonian ((2|) in the single-excitation 
subspace thus finding its eigenvalues Ej and the corre- 
sponding normalized eigenstates \Ej). We discuss the 
system dynamics for different energy regimes and con- 
sider no detuning, A = 0, between photonic and atomic 
modes. We are mostly interested in the analysis of pho- 
tonic and atomic excitation occupation probabilities in 
different cavities for a given initial state, namely \ip(0)), 
through the time. The state at time t is given by 
\${t)) = U(t)\Tp(0)), where U(t) = e~ im is the quan- 
tum time evolution operator. If the initial state is pre- 
pared in a equal superposition of photonic and atomic 



A. Large hopping regime 

First, we investigate the large hopping regime (k ^> 0). 
In this case, the atom-cavity interaction slowly takes 
place so that the photonic excitation propagates practi- 
cally free, while the atomic one does not "move" in short 
time scales. Remember that only the photons can hop 
while atoms remain stationary. The atomic component 
dynamics is then correlated to the photonic one, although 
not exactly the same, except when the atom-cavity cou- 
pling /3 dominates the hopping strength as we will see 
later. 

The AN has a 27r/3 rotational symmetry, hence there 
are several subsets of equivalent nodes (see Fig. [2]) like 
{1,2,3}, {5,6,7}, {8,10,11,13,14,16}, and so forth, in 
which |^fe(0)) would lead to the same dynamics. Like- 
wise, for the initial state set at the hub node (k = 4), 
the occupation probabilities ^^(t) are equal for each 
subset. 

Figs. [3] and 2] show the time evolution of the pho- 
tonic occupation probability 7r? fe (t), for different initial 
single-cavity eigenstates \ipk(0)), in the n — 3 and n = 4 
AN, respectively. Although exhibiting a rather complex 
dynamics, we can identify some general aspects of the 
propagation. The most remarkable feature is periodic 
transition between localized and extended states. The 
return probability 7rjy. (t) can reach the highest values so 
that the excitation is most likely to be found in the ini- 
tial node. This localization is stronger for k = 4, the hub 
node. Also, it seems that 7r^(t) oscillates faster for the 
initial state set at nodes with a higher degree of freedom 
7. In Fig. [5]we plot 7rf 4 (t) (solid line) and also the atomic 
component 7rf|(t) (dotted line) time evolution in detail. 
As we said before the atomic mode barely moves, ex- 
cept for small oscillations (see inset of Fig. [5]) due to the 
photonic component behaviour. Since the atom-cavity 
interaction process has much lower associated energies, 
only the cumulative effect of these interactions in longer 
time scales may significantly vary 7r|^ (t ) . Therefore, at 
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FIG. 3. (Color online) Time evolution of the photonic exci- 
tation n^(t) through the n = 3 AN (16 nodes) for the large 
hopping regime. The vertical axis represents the occupation 
probability of a particular node (cavity). The single-cavity 
initial state \tp k (0)) = \k) CS> (\g, 1) + [e, 0})/y/2 is prepared at 
nodes k = 3, 4, 5, and 15. System parameters are k//3 — 10 3 
and A = 0. The node numbering is available in Fig. [2] 



the large hopping regime we shall focus on the photonic 
mode propagation dynamics thus aiming to explain the 
quantum transport properties of the AN. 

Now let us turn our attention to the probability of 
finding the excitation in nodes other than the initial one. 
In general, such probabilities are very small compared 
to TTkk(t), although they can reach relevant values in 
some specific nodes. The overall probability distribution 
strongly depends on the AN generation and the initial 
conditions. In Fig. [31 by examining irik{t) for fc = 3, we 
see that a significant amount of probability flows to the 
other nodes belonging to the same degree group, that is 
1 and 2 (check Fig. [5] for node numbering). The same 
applies to 7r 65 (i) and ^75 (t) when k = 5. For k — 15, the 
nodes 10 and 11 even belonging to the same degree group 
are barely populated. It suggests that sharing the same 
degree does not necessarily imply in privileged network- 
ing. It turns clearer in Fig. [4] for k = 18 and 21, where 
the excitation is mostly spread to a few nodes within the 
subnetwork [32| generated by the new corner nodes 1,2, 
and 4. When k = 5 for n ~ 4 (Fig. @|, the excitation still 
roams between nodes 6 and 7, but now avoiding trans- 
port to higher generation nodes. For k = 15 in Fig. 21 
there is a relevant probability of finding the excitation 
at nodes 20, 21, and 22. These nodes belong to the sub- 
network formed by nodes 1, 2, and 5, where the node 15 
play the role of hub node. These observations indicate 
that these subnetworks has a fundamental part on the 
system transport dynamics for higher AN generations. 



FIG. 4. (Color online) Same as Fig. [3] but for n = 4 (43 
nodes) and k = 5, 15, 18, and 21. 
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FIG. 5. (Color online) Photonic (solid line) and atomic (dot- 
ted line) return probabilities for the initial state prepared in 
|^/)4(0)) at the n — 4 AN. The inset shows small oscillations 
in 7i~44(f) induced by the photonic mode propagation. System 
parameters are k//3 = 10 3 and A = 0. 



A better way to visualize how the excitation probabil- 
ity distributes between the nodes for given initial con- 
ditions is by calculating the long time average of ir^ (t) 
defined as 



<u(t)dt 



= Y J KE l ^E. ] ){t. / g\E l ) (^|Vfc(0)) 



(7) 
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FIG. 6. (Color online) Photonic long-time average occupation probabilities for (a) n = 4 (43 nodes), (b) n = 5 (124 nodes), 
and (c) n = 6 (367 nodes). System parameters are k//3 = 10 3 and A = 0. The plot shape is symmetrical since x^k — x£ 
for every k and I, and the highest values are i.e., the average return probability. Although we have not provided the 
node numbering for nodes from generations higher than n = 4, the purpose here is to give an idea of how the subnetworks 
and nodes from different generations are connected among themselves. In each generation above, the darkest branch represents 
the average excitation transport probability between nodes belonging to the current generation and to the previous one. The 
square boundaries point out connections between nodes having the same degree. 



where 5(J5j - Ej) = 1 for E. t = Ej and 6(Ei - Ej) = 
else. Fig. [5] shows the distribution of this average for dif- 
ferent AN generations. First of all, notice that Xtk — x\t 
for every k and £ (since the evolution is unitary), hence 
this quantity provides us a clearer view on the way the 
nodes stablish communication channels with each other. 
As it must be, the higher averages are Xfefc> however, by 
looking at the other connections we can identify several 
zones (square boundaries) and branches apart which, in- 
terestingly, present self-similarity as n increases. The 
branches connect generations with each other while the 
zones connect nodes from the same generation, i.e., with 
same degree 7. The strongest (darkest) branch is always 
the one connecting nodes from the current generation to 
the previous one. It mostly represents the links between 
the corresponding hub nodes from the 3™ -2 subnetworks 
having 7 nodes, to its three nearest-neighbour nodes with 
7 = 3. As n increases, the previous branches still remain 
but with lower intensities since the aforementioned hub 
nodes now have a higher 7. Also, nodes from "distant" 
generations seem to be barely connected. This is because, 
when preparing the initial state in nodes from a current 
generation, the excitation is mostly roams between the 
initial node, other nodes with same degree belonging to 
local subnetworks, and its corresponding hub node from 
the previous generation. Notice in Fig. [5] that connec- 
tions between nodes from the same generation are only 
relevant within subnetworks composed up to 16 nodes. 
For instance, at a given generation n there are 3™~ m 
subnetworks having (3 m + 5)/2 nodes. Needless to say, 
these degree interconnections also become weaker as gen- 
eration increases. 

Before discuss other energy regimes it would be use- 
ful to explain the reason for why the state always carries 
relevant information about its initial conditions by eval- 



uating the participation ratio of eigenstates defined as 



SJ N(n) ' v ; 

E E IM£;>I 4 

s={g,e\ 1=1 

for a given AN generation n. This quantity can assume 
values within 1, for completely localized states, and D, 
for extended states ((i,s\Ej) = l/V~D for all i). Fig. [7] 
shows the relationship between £j/N and the probability 

of finding the system in each eigenstate, \(Ej\ipk(0))\ 2 , 
for a given initial condition \tp k (Q)) [Fig. [Tfa)], and the 
associated eigenvalues Ej [Fig. E[b)], for the n = 4 AN. 
When we prepare [^(0)) in the hub node k = 4 (with 
degree 7 = 24), the two most localized eigenstates (with 
£,j/D « 0.05), one being related to atomic excitation 
|4, e) (Ej/ k sa 0) and the other to the photonic excitation 
|4, g) (the highest eigenvalue), become more populated 
thus strongly leading the system dynamics. Only a few 
other eigenstates, not much localized, gives a relevant 
contribution to the dynamics. Since we are considering 
the large-hopping regime, the atom-cavity interaction en- 
ergies are very low so we clearly identify the middle range 
of eigenvalues ss [see Fig. [Tfb)] as being those belong- 
ing mostly to the atomic component of the system. In 
Fig. [T](a) , as we set ^(O)) for nodes having lower de- 
grees of freedom, for instance, k = 3, 5, and 21, where 
7 = 17, 12, and 3, respectively, the contribution from 
localized eigenstates decreases and then more extended 
states start to participate on the dynamics. Also, the 
eigenstates which mostly contribute to 7r^(t) has lower 
associated energy levels for smaller 7. It shows that the 
oscillation period of the return probability depends on 7, 
as we have already mentioned. 

More details about localized and extended states on 
the AN and its relationship with the structure charac- 
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FIG. 7. (Color online) (a) The normalized participation rate 
of eigenstates £j/D (solid line) and the probability of find- 
ing the system in eigenstate \Ej), that is \(Ej\ipk(0))\ , for 
k = 4, 3, 5, and 21. System parameters are n = 4 (43 nodes), 
k//3 = 10 3 , and A = 0. The eigenstates are labeled accord- 
ing to increasing energy values and, within each degenerate 
group, by increasing values of £j\ (b) The corresponding en- 
ergy values (in units of k _1 ). The energy levels around zero 
is shown on the inset. Actually, the ones lying at the middle 
range (sa 0) are not all degenerate. These energy levels arise 
from the atomic modes and since s » (1, its energy scale is 
very small compared to the photonic counterpart. 

teristic spectrum and symmetry properties can be found 
in Ref. [281 ]. where they considered a simple tight-biding 
electronic model. The tendency of the excitation for stay- 
ing in the initial node has been found for the continuous- 
time quantum walk in the same structure [3fj| . and also 
in a small- world network defined by adding random links 
in a one-dimensional ring of nodes [241 ] . By increasing 
the number of additional bonds, they have shown that 
this effect becomes stronger thus avoiding dispersion. 

B. The JCH and strong atom-cavity coupling 
regimes 

In what follows we consider the situation where the 
coupling between cavities and the atom-cavity interac- 
tion strength are of the the same order (k « 0). That is 
the energy regime where the JCH model displays its full 



nature. In this the photonic excitation propa- 

gates through the network it exchanges energy with the 
atomic modes in such a way that the whole system starts 
to roam between photonic and atomic degrees of freedom. 
Each cavity provides its contribution to mostly either one 
of these two components, depending on the system ini- 
tial conditions and the characteristic single-cavity energy 
levels in a rather complex way. Since a fully localized ini- 
tial state of the form |?/>fc(0)) mostly leads the excitation 
to return to the initial node, we presume that the oscil- 
lation rate of ituu pe will be primarily relevant for driving 
the whole system among the atomic and photonic degrees 
of freedom. Fig. [8] gives an idea of how the initial state 
may significantly changes the way both modes contribute 
to the system dynamics in such energy regime. In Fig. 
[5£a) we show the long-time average occupation probabil- 
ity distribution Xik & wnen IV'fc(O)) is set at the hub node, 
k = 4. Notice that xte can ^ e evaluated by the same way 
as Eq. [7]by projecting the eigenstates on \£, e), instead of 
\£, g). We can see that some groups of nodes have higher 
X? 4 values (including the hub node), others mostly con- 
tribute to xfli an d a fe w practically contributes with the 
same amount for both components. In Fig. [5Jb) we plot 
the time evolution of the probability of finding the pho- 
tonic and atomic excitation at any cavity by summing 
7r|| pe over £. In this case, the system notoriously fluctu- 
ate between both components. It results mainly from the 
fast oscillation rate of 7144 (since the hub node has the 
higher 7) combined with the strong localization proper- 
ties. The average values are X44 ~ 0.26 and Xm ~ 0.14 
hence representing a great difference. Figs. [SJc) and|U[d) 
show the same quantities but now for k = 18 (7 = 3). 
The higher averages x^g 6 does not differ much for both 
types, thus resulting in small variations of n eis e ■ 

In summary, at the large hopping regime, k» /?, the 
interplay between atomic and photonic degrees of free- 
dom is minimal so that the atomic excitation only signif- 
icantly propagates in larger time scales. Thus the cavity- 
cavity coupling parameter, k, allows the control of the 
atomic mode velocity. When the JCH regime (k ~ 0) 
takes place, atomic and photonic degrees of freedom start 
to interfere with each other, thus resulting in different dy- 
namics . We can keep on decreasing k//3 until reaching 
the strong atom-cavity coupling regime, where photonic 
and atomic components have exactly the same dynam- 
ics. This regime transition is clearly shown on Fig. [S] 
where we plot the time evolution of the photonic [Fig. 
[5Ja)] and atomic [Fig. Wft>)] return probabilities for 
a large range of k//3 values from 10 3 to 10~ 3 . At the 
k ft regime, the dynamics for both atomic and pho- 
tonic modes behaves exactly like at the large hopping 
regime but now creating new single-cavity eigenstates (if 
the initial state is prepared in \ipk(0))) as passing by other 
cavities, and propagating at half the speed of the pho- 
tonic excitation at the k ^> ft regime [18| . Due to the 
fact that the atoms inside the cavities does not actually 
move, there is a small delay, namely 8, which tends to dis- 
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FIG. 8. (Color online) Figures (a) and (c) shows the photonic and atomic long-time average occupation probability distribution 
for k = 4 and k = 18, respectively, and (b) and (d), the corresponding time-dependent probability of finding both excitation 
types at any cavity, £\ 7r ffe Pe - System parameters are n = 4, k//3 — 1, and A = 0. The return average probabilities are 
X p 4 1 » 0.26, X 44 « 0.14, x?8i 8 « 0.103, and X iLs » 0.105. 



turb the single-cavity eigenstates, but S — > as k//3 —> 0. 
By the same way, at the large hopping regime, 6 — > oo as 
k//3 — > oo thus freezing the atomic component. 



IV. CONCLUSIONS 

Even at the single-excitation subspace, the JCH model 
provides a rich dynamics. In addition, the Apollonian 
network has ubiquitous transport features. We studied 
the propagation dynamics by considering the initial state 
as an even superposition of atomic and photonic modes 
fully localized in some specific node. At the large hopping 
regime, where practically only the photonic component 
moves, we focused the analysis on the quantum transport 
dynamics in the AN. In such structure, the excitation is 
most likely to be found at the initial node due to the 
characteristic network spectrum properties. By evalu- 
ating the long time average photonic occupation prob- 
ability for different AN generations, we have discussed 
the interplay between the nodes and subnetworks thus 
showing their role on the system dynamics. At a given 
generation n, the nodes with 7 = 3 are mostly connected 
with its correspondent hub node from the previous gener- 
ation and with some other nodes having the same degree 



within local subnetworks. These properties are also valid 
for any analogue tight-binding model in the same AN. 



At the JCH regime, where the hopping term and the 
atom-cavity coupling strength are of the same order, we 
have shown that each cavity gives its own contribution for 
the atomic and photonic modes thus inducing the system 
to roam between these two degrees of freedom. When 
the atom-cavity coupling term dominates, the excitation 
propagates maintaining the even superposition of atomic 
and photonic modes (the single-cavity eigenstate). 



Even though high-fidelity excitation transport pro- 
tocols in an one-dimensional chain of coupled cavities 
has been proposed [181 Il9f , such regular structures have 
strong dispersive properties, in contrast to the AN, where 
we have shown that the excitation actually avoids prop- 
agating to some specific nodes, even when it is a first- 
neighbour node, without any coding-decoding protocol. 
In fact, we have only dealt with the most simple fully lo- 
calized entangled initial state. Obviously, there are sev- 
eral possible combinations of initial states that should 
generate different and much more complex dynamics. 
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(a) (b) 




lOg 10 (K/P) lOg 10 (K/P) 

FIG. 9. (Color online) Comparison between atomic and photonic propagation dynamics from the large hopping regime k//3 = 
10 3 to the strong atom-cavity coupling regime n/fi = 10~ 3 . In (a), the photonic return probability to the hub node, 71-44 > an d 
in (b) the atomic counterpart. System parameters are n = 4 and A = 0. The middle range of ft//3 values is where the JCH 
model plays its main role. 
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